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ABSTRACT 

With gravity, ionization, and radiation being considered, we perform 2.5D 
compressible resistive MHD simulations of chromospheric magnetic reconnection 
using the CIP-MOCCT scheme. The temperature distribution of the quiet-Sun 
atmospheric model VALC and the helium abundance (10%) are adopted. Our 
2.5D MHD simulation reproduces qualitatively the temperature enhancement 
observed in chromospheric microflares. The temperature enhancement AT is 
demonstrated to be sensitive to the background magnetic field, whereas the 
total evolution time At is sensitive to the magnitude of the anomalous resis- 
tivity. Moveover, we found a scaling law, which is described as AT / At ~ 
hh^^'^ B'^'^rjo^-^^ . Our results also indicate that the velocity of the upward jet 
is much greater than that of the downward jet and the X-point may move up or 
down. 

Subject headings: MHD — Sun: activity — Sun: chromosphere — Sun: flares 



INTRODUCTION 



Magnetic reconnection plays a very important role in solar flares, corona mass ejections. 



and other solar activities. During the solar minimum. 



on the solar small-scale activities su ch as microflares flOiu et al. 



Nine 



2008 



Brosius fc Holman 



many authors 



20091). EUerman b ombs ( iFang et al. 



20081 ). chromospheric jets (INishizuka et al.ll2008l ). and so on. 



2004 



bcus their attention 



Fang et al. 



3 



2006a: 



2006bl: IWatanabe et al. 



Microflares, or subfiares, or bright points, are small-scale and short-lifetime solar 
activities. The size of microflares ranges from several arcsecs to about 20 arcsecs, the 
duration and the total energy can be 10 — 30 minutes and 10^^ — 10^^ ergs, respectively 
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( IShimizu et al.ll2002l : iFang et al.ll2006al ). The most distinctive feature in their visible spectra 



is the faint emission at the center and the brig. 



hnes such as 



ia hne. Soft X-ra y ( 



1984 ). EUV flPorter et al. 



1984 



itening at t he wing of some chr omospheri c 



Golub et al 



1974 



Emslie fc Noyes 



9771 ). hard X-ray 



(HXR 



19781 ) and microwave (IGary fc Zirin 



Lin et al 



19881 ) emissions have also been observed in some microfiares. It is an interesting fact that 



there are s ome c orrelations among the emissions at different wavelengths. For instance, 



Qiu et al 



(120041 ) found that about 40% of microfi ares show H XR emissions at >10 keV 



and microwave emissions at ~10 GHz. Recently, 



Ningl (|2008[ ) found that there is a 



correlation bet ween the power-law index of the HXR spectrum and the emission measure 



of microfiares. iBrosius fc HolmanI (j2009l ) found that the studied microfiare is bright first 
in the chromospheric and transition region spectral lines rather than in the corona, which 
is consistent with the chromospheric heating by nonthermal electron beams. All these 
evidence indicates that microfiares are related to the nonthermal processes driven by 
magnetic reconnection. 

Some microfiares might result from the magnetic reconnection in the corona. However, 



reconnection in the chromosphere (iLiu et al. 



LKKestmg \ 
(Liu et al. 


nax s 
2004; 


ome oiner 
Xia et al. 


micr( 
2007) 



Tang et al. 



( I2OOOI ) found that in many cases emerging fiux occurred about 5 — 30 minutes before the 



microfiares, where the emerging fiux may c ollide with the pre-existing magnetic field at the 



chromospheric height. 



Brosius fc HolmanI ( 120091 ) also stated that magnetic reconnection 



in the chromosphere could be a plausible mechanism for triggering the chromospheric 
microfiare. 



Chen et al. 



(120011 ) made 2.5-dimensional (2.5D) numerical simulations of chromospheric 
magnetic reconnection in order to study Ellerman bombs and type II white-light flares. 
However, in the simulations, they made some simplifications, such as the omissionof 
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the gravity. In this paper, in order to improve their results, we perform 2.5D 
magnetohydrodynamic (MHD) simulations using the CIP-MOCCT scheme, with gravity 
being included, and try to make some comparisons with the observations of chromospheric 
microflares. In the next section, the numerical method is given in detail. Our results, which 
include the dynamic process, parameter dependence, a scaling law and a comparison with 
a semi-empirical model, are described in Section [31 Discussion and summary are given in 
Section HI 



In our simulation, we solve 2.5D compressible resistive MHD equations with ionization, 
radiation, and uniform gravitational field being considered. There is no thermal conduction 
in our simulation since thermal conductivity is very small in the chromosphere. The MHD 
equations are given as follows: 



2. 



NUMERICAL METHOD 



2.1. Basic Equations 



dp 
'dt 



+ V ■ (pv) = , 







(3) 



de 
di 



+ V ■ [v(e + 



+ p)- — B(B ■ v) + (^V xB)xB]=pg-Y + H-Rr , (4) 



2po 
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where eight independent variables are the density (p), velocity {vx, Vy, Vz), magnetic 
field {Bj., By, Bz), and temperature (T). The thermal and ionization terms are included 
into the total energy density e = p/(r — 1) + + B'^/2fiQ + XH^e, where xh is the 

ionization potential of hydrogen atoms and /iq is the vacuum permeability. Considering 
the helium abundance (10%), we take p = lAmnnH, where uh is the number density 
of hydrogen. With partial ionization considered, the total gas pressure is expressed as 
p = [IAuh + ne)kBT, where kB is the Boltzmann constant and the number density 
of electrons, rig, is de duced from the mo dified Saha and Boltzmann formula for a pure 



hydrogen atmosphere (iGan &: Fang 



1990|): 



He = (V0' + 4nj^0-0)/2 



(5) 



1 ,27rme/tRT, -^h 



(6) 



2T 

hi = — e^-RT^T 
Tr 



(7) 



19731 ). In Equation 



where ho is the Planck constant. We take Tr = 6000 K (iBrown 
(jl]), Rr and H are the radiative loss and the pre- heating term, respectively. Radiative 
loss is important in the chromosphere and the photosphere. Generally speaking, it 
should be calculated by solving simultaneously the statistical equihbrium and radiative 
transfer equations along with the MHD equations. However, it is not feasible to deal 
with it in our 2.5D numerical simulation. Instead, we fit the non-local thermodynamic 
equilibrium (non-LTE ) radiative losses in the VALC model for the quiet- S un atmosphere 



fIVernazza et al. 



198ll ) and the Fl semi-empirical mod el for the weak fiare (IMachado et al. 



19801 ). and then modify the empirical formula given by 



Gan fc Fang] (jl990l ) in order to make 



- 6 - 



it more suitable for the small scale activities such as microflares: 

Rr = nHnea{y)f{T) , (8) 

a{y) = 10(1-745x10-3,-4.739) ^ g_Q X ^Q-2g-3.701xl0-, ^ ^g) 

/(T) = 4.533 X 10-23(^)2-«^^ , (10) 

where a{y) and /(T) are functions of y (the height from tsooo = 1 in the photosphere) and 
the temperature T, respectively. The unit of y is km. The suitable range of the height 
and the temperature in Equation ([8]) is < 2000km and T < 10^ K, respectively. The 
pre-heating term is given by if = nH{neOif)t=o-, which balances the radiative loss before 
magnetic reconnection. 

2.2. Normalization and Parameters 

In order to non-dimensionalize the MHD Equations ([I]) - (HI) in Section 12.11 the 
variables are normalized by the quantities given in Table [H In this table, mj;^ is the 
proton mass. The dimensionless MHD Equations, which have a standard form of advection 
and non-advection terms, are numerically solved with the CIP-MOCCT scheme (Kudoh 
1999). Besides the acoustic transit time scale to used in the nondimensionalization, another 
practical time scale is the Alfven transit time ta = Lq/va, where va is the Alfven velocity 
at the height of the reconnection X-point. In our simulations, two cases are considered, 
where the magnetic reconnection X-point is set at two different heights, i.e., 500 km and 
1000 km. The corresponding Alfven velocities are 18.2 km and 11.2 km s^^, while the 
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corresponding Alfven transit time is 109.9 s and 177.6 s, respectively. 



2.3. Initial and Boundary Conditions 



The domain of the numerical simulation is -1000 km < x < 1000 km, < y < 2000 
km, as shown in Figure [1] The initial magnetic configura tion is a force-free field with a 
current sheet located along the y-axis, which was used by 
as follows: 



Chen et al. 



( 119991 ) and rewritten 



B. = 0, 



:iii 



B„ 



R a; 



\x\ > Ah 



i?osin 2^ < Ah 



(12) 



|x| > Ah 

i?ocos^ |x| < Ah 



(13) 



where Bq and Ah are the background magnetic field and the half width of current sheet, 
respectively. In order to make our simulation more realistic, we use the temperature 
distribution of the quiet-Sun atmospheric model VALC as shown in Figure [2] (left panel). 
Using the hydrostatic equilibrium equation and the modified Saha and Boltzmann formula 
for calculating the number density of electrons iteratively, we get the initial number density 
of hydrogen (uh) and electrons (ue) as shown in the middle and the right panels of Figure 
El respectively. From these figures we can see that our initial model is very close to the 
VALC mo del. In order to initiate the magnetic reconnection, we assume that an anomalous 



resistivity (j Cramer fc Donnelly 



19791), with the form rj = tjq cos[x7r/2A/i] cos[{y — hr)TT/2Ah], 



is localized in a small region |a;| < Ah, \y — hr\ < Ah (the small rectangle in Figures (H [7| 
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and [T2|) . where hr stands for the height of the reconnection point (X-point) in the y-axis, 
Ah is the half width of the resistivity region and rjo is the resistive amphtude. It is noted 
that the half width of the resistivity region is taken to be the same as the half width of the 
current sheet. 

Owing to the symmetry, the calculation is performed only in the right half region, with 
the symmetry boundary condition being set on the left-hand side. The right and the upper 
boundaries are treated as open boundaries. Line-tying conditions are applied to the bottom 
boundary. It is noted that p should satisfy the hydrostatic equilibrium {dp/dy + pg = 0) 
to reach a balance between the gravity and the pressure gradient force at the upper and 
the lower boundaries. The calculation domain is discretized into 800 x 400 grid meshes. 
The grid sizes are Ax = 1.25 km and Ay = 5 km, which are much less than the pressure 
scale height (> 100 km). With gravity, ionization and radiation being considered, and 
with heat conduction being neglected, the numerical simulations are performed with the 
CIP-MOCCT scheme. 



NUMERICAL RESULTS 



3.1. Dynamic Process 



Generally speaking, the height of the X-point 



( iTakeuchi &: Shibata 



2001 



von Rekowski &: Hood 



may change as the reconnection proceeds 



20081 ). In our simulation, we found that 



the X-point could move up or down by 10 — 40 km, but still within the resistivity region. 
For this reason, we do not need to change the position of the resistivity region. In order to 
illustrate the dynamic process, this subsection describes the simulation results in a typical 
case, where the background magnetic field Bq is taken to be 50 G, the anomalous resistivity 
?7o be 17.7 Q m, the X-point height be 1000 km, and the half width of the resistivity 
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region Ah be 100 km. Figure [3] depicts the results of the simulation at different times. In 
this figure the color stands for the temperature, solid lines for magnetic field and arrows for 
velocity. 

As the anom alous resistivity sets in, we have a fast reconnection process similar to 



that predicted by iPetscheld (jl964j ). Figure H] shows the evolution of the reconnection rate, 
which is defined as \vin\/vA, where Vm is the maximum inflow velocity along the x-direction 
at the height of the X-point. It is found that the simulation stops at t = 1.35r4 (148.2 s), 
probably due to some numerical instability. The reconnection rate reaches the maximum, 
0.11, at t = I.Ita (120.8 s). After that, the reconnection rate remains almost constant 
around 0.1. The lifetime of the reconnection simulated in this case (148.2 s) is shorter 
than the typical lifetime of microfiares, i.e., 10 — 30 minutes. Maybe our results can only 
represent the pre-fiare phase and the rise phase. 

During the reconnection, the temperature enhancement in the resistivity region is 
about 2000 — 3000 K as shown in Figure O Accelerated by the magnetic tension force, 
the hot plasma together with the frozen-in reconnected field lines are ejected from the 
resistivity region to form upward and downward reconnection jets. The maximum velocity 



-1 



of the upward and downward jets are about V^p = 23.0 km s~^ and Vdown = 9.4 km s 
respectively as shown in Figure El Because of the sharp density gradient in the lower 
atmosphere, the speed of the upward jet is faster than that of the downward one. Compared 
with the local Alfven speed {Va = 18.2 km s^^) at the height of 1000 km, we find that 
the velocity of the outflow is approximately the local Alfven speed in the inflow region. 
The downward jet can only reach the height of 600 km, after which it is slowed down 
significantly. Furthermore, the reconnection X-point also rises by about 40 km during the 
reconnection process. 

The hot plasma jet ejected along the y-axis is mainly heated by the Joule dissipation 
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in the resistivity region. However, we can not ignore the heating by the slow-mode shocks. 
In our simulation, the heating occurs not only in the resistivity region but also outside the 
region as shown in Figure O which shows the current density (z-component) distribution in 
the left panel and the temperature distribution in the right panel. It can be seen that a pair 
of slow-mode shocks, which are characterized by the high current density, are formed at the 
interface between the reconnection inflow and outflow, across which the plasma temperature 
increases drastically. It indicates that besides the Joule dissipation, the slow-mode shock 
may be another effective way to heat plasma. 



3.2. Parameter Dependence 

In our simulation, the most important parameters are the background magnetic field, 
the anomalous resistivity, and the height of the X-point (i.e., the reconnection site). In 
order to show the influence of different parameters on the results, we performed extensive 
simulations by changing one parameter at a time with others being fixed. 

Figure [8] shows the evolution of the magnetic reconnection rate for different initial 
background magnetic fields and anomalous resistivities, with the height of the initial 
X-point being 1000 km. It can be seen that the magnetic reconnection rate can reach the 
maximum much faster when we set a stronger initial magnetic field. It is also found that 
the large anomalous resistivity can shorten the evolution time from the initial stage to 
the time when the reconnection rate reaches the maximum. It should be noted here that, 
although the resistivity in the chromosphere is ~ 3 orders of magnitude larger than that 
in the corona, it is still too small from the simulation point of view. A so small resistivity 
would lead to a very thin diffusion layer near the reconnection X-point, which could not be 
resolved by the current numerical mesh. Therefore, the values of the anomalous resistivity 
we used are several orders of magnitude larger than the classical Spitzer resistivity (e.g.. 



Kovitya & Cram 1983). 



Figure IH] shows the temperature distributions at the peak time of the magnetic 
reconnection rate. Our results indicate that the temperature enhancement is sensitive to 



the magnetic field. However, in all cases our results are consistent 



qualitative 



y wit 



1 the 



2006ah . A 



observed temperature enhancement in chromospheric microfiares (iFang et al. 
quantitative comparison between the simulation and the observation will be described in 
Section 3.4. 



3.3. Scaling law 

The ranges of the variables such as the background magnetic field (Bq), anomalous 
resistivity {rjo) and the number density of hydrogen {nn) are very wide in observations. 
As a result, the maximum temperature enhancement (AT) and the reconnection lifetime 
(At) may change greatly. In order to find a scaling law to relate AT and At to the 
above-mentioned 3 variables, we change the variables one by one and perform a series 
simulations. As illustrated by Figure [101 the dataset can be well fitted by the following 
formula: 



From the energy point of view, as the reconnection proceeds, the Joule dissipation is 
consumed partly by the plasma radiation and ionization, with the rest being converted into 
the plasma heating. Therefore, we have 



nnkAT = — At - nHa{y)nJ{T)At + nHa{y)[nef{T)]t=QAt - [n^ - {n^)t=o]XH , (15) 
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where the expressions of rig, a{y), f{T) are given by Equations ([5]), ([9]), and f lTOj) . The 
items on the right-hand side of the Equation f[T51) are the Joule dissipation, the radiative 
losses, the pre-heating and the ionization energy, respectively. If we assume that the plasma 
radiation and ionization are negligible. Equation (fT5|) can be rewritten in a simple form as 



nnkAT ^ -At = ^^^^^^At ~ ^-^At , (16) 

where /zq is the vacuum permeability, L is the scale of the reconnection region. Since in our 
simulation L is fixed as a constant, we can get a simple theoretical formula as follow: 



At uh 

It can be seen that the formula given by our simulations is similar to the formula 
f lT7|) except the relatively bigger deviation of the power index of rtu. The deviation is 
probably because the radiative loss and ionization items, which are related to hh but not to 
Bq and ?7o, have been considered in our simulations but are neglected in the analytical form. 
A further study will be done in order to find a more realistic scaling law in the future. 



3.4. Comparison with a Bright Microfiare Event 

Fang et al. (2006a) analyzed five well-observed microfiares and computed the thermal 
semi-empirical models for two typical microfiares, one was bright and the other was faint. 
In oder to get the semi-empirical models, they first assume the temperature and turbulent 
velocity distributions, and then derive the density distribution. Based on these distributions 
they can calculate the Ha and Ca II 8542 A line profiles. The atmospheric model is adjusted 
until the calculated line profiles to match the observed ones. The temperature distribution 
of the bright microfiares model is shown in the middle panel of Figure [TTl In order to 
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reproduce the temperature enhancement in the model, we set the initial background 
magnetic filed to be 300 G, the anomalous resistivity to be 159.3 fl m, and the initial 
reconnection height to be 500 km. Since the magnetic field is strong, the total evolution 
time is short. The evolution of the reconnection rate, the temperature enhancement, and 
the velocity along the y-axis {vy) are shown in Figure [Til 

After t = 0.2ta, the reconnection rate reaches about 0.1. We choose the temperature 
distribution at 0.24 ta to compare with the semi-empirical model. It is shown that there 
are two peaks in our simulation results. The smaller one is located at the height of 350 
km and the bigger one is almost the same as the semi-empirical model and located at the 
height of 600 km. The possible reason for this phenomenon is slow-mode MHD shock 
heating as shown in Figure [121 Below 800 km, the temperature profile of our simulation is 
close to that of the semi-empirical model. However, above 800 km, the temperature in the 
semi-empirical model becomes much higher than that of our simulation. The discrepancy 
might result from the fact that our simulations are performed for a relatively short period 
such that the reconnection jet has not yet reached the corona. Therefore, our result can 
only give a reasonable explanation for the temperature bump of the microfiare mainly in 
the rise phase. 

The right panel of Figure [TT1 shows the one-dimensional plot of Vy. From this figure it 
can be seen that the height of the reconnection X-point, which is characterized by the null 
velocity, has dropped from 500 km at t = to 450 km at t = 0.24r^. The velocity of the 
upward fiow is about 11.1 km s^^, which is much larger than that of the downward fiow 
(1.8 km s~^). The downward jets can only reach the height of about 200 km. 
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DISCUSSION AND SUMMARY 



More and more theoretical works have indicated that r nagnetic reconnection in the so lar 



lower atmosphere can produce chromos pheric microflares (iTandberg-Hanssen fc Emslie 



1988 



Liu et al. 



bombs ( 



2004 



Fang et al. 



Fang et al 



2006b 



Ding et al. 



2006a), as well as other small activities like EUerman 



1998 



Watanabe et al. 



20081 ). With gravity, ionization 



and radiation being considered, we performed 2.5D MHD simulations. Si nce we have used 



the t emperature distribution of the quiet-Sun atmospheric model VALC (IVernazza et al. 



198ll ) and considered the helium abundance (10%), our simulations are realistic to some 



extent. 



In our simulation, we calculate the magnetic reconnec tion rate by using the de finition 



R = Vin/vA, rathe r than the formula R = dtlj /dt used by 



Chen et al. 



R = IV'Jzl used by 



(1999 



20011 ) and 



Yokoyama fc Shibatal (120011 ). where ip means the magnetic flux function 



and Jz stands for the 2;-component of current density. If we use the formula R = dip/dt 
or R = {fjJzl, the reconnection rate has a large pulse, which can b e greater than 1.0 at 



the be ginning of the evolution. It is non-physical, as mentioned by 



Yokoyama &: Shibata 



( 120011 ). However, it should be noted that all of these methods can give a similar result after 



enough time of evolution. 



In the subsections 13.11 and 13. 4[ we studied the magnetic reconnection with the X-point 
at different heights, e.g., 500 km and 1000 km. Our results, in the cas e of 500 km, can w ell 



2006al ). 



reproduce the temperature enhancement in the semi-empirical model (iFang et al.l 
At the heights of 500 km and 1000 km , the temperature enhancement is 1000-2000 K and 
2000-3000 K, respectively. In the case of 1000 km, we found that the height of X-point can 
move up by tens of kilometers. This is different from the results of the case of 500 km, 
where the X-point moves down. It is noted that the temperature enhancement we got in 



our simulations is several thousand Kelvin, and our results can only represent the rise phase 



- 15 - 



of chromospheric microflares since our simulations were performed only until the upward 
reconnection jet reaches the bottom of the transition region. Of course, those microflares 
with strong EUV, X-ray, and microwave emissions may occur at the lower corona. In this 
case, the emissions in X-rays, EUV and microwave can be naturally explained. 

It is worth noting that our simulations have some limitations. Even though we limited 
the computational height below 2000 km, the plasma /3 at the bottom boundary is still 
six orders of magnitude larger than that at the top boundary. It causes some numerical 
instability in our simulations. Besides, our initial magnetic configuration is relatively 
simple. Moreover, as it is well known, the current sheet should be very thin in the real 
situation, e.g., less than hundreds of meters, while in our simulation the minimum grid sizes 
are Ax = 1.25 km and Ay = 5 km, which are too large to simulate the current sheet. All 
these contribute to the fact that our simulations can not reproduce the real observations in 
details. Further improvement of this work is expected. 

In summary, we give the conclusions as follows: 

1. Our 2.5D MHD simulations can reproduce the temperature enhancement in 
chromospheric microflares qualitatively. The temperature increase in the cases when the 
reconnection point is at the height of 500 km and 1000 km can reach 1000-2000 K and 
2000-3000 K, respectively. 

2. The free parameters in our 2.5D simulation are the background magnetic field (Sq) 
and the anomalous resistivity (r^o)- We have found that the temperature enhancement is 
sensitive to the background magnetic field, while the total evolution time is sensitive to the 
magnitude of anomalous resistivity. 

3. Our simulation results indicate that the velocity of the upward jet is much larger 
than that of the downward jet, and the X-point may exhibit downward or upward motions. 



4. We have performed a parameter survey, and found that the temperature 
enhancement and the total evolution time in the chromospheric reconnection follow the 
scaling law: 



AT 



3.7( 




Bo 



) 



2.1 




-1 



(18) 



At 



3.4 X 10^9 m-' 



25 G 



17.7 n m 



The computations were done by using the HP ProLiant BL260C G5 Blade system 
at the Department of Astronomy of Nanjing University of China. This work is supported 
by the National Natural Science Foundation of China (NSFC) under the grant numbers 
10221001, 10878002, 10403003, 10620150099, 10610099, 10933003 and 10673004, as well as 
the grant from the 973 project 2006CB806302. 
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Table 1: Normalization Units 



Variable 


Quantity 


Unit 


Value 


x,y 


Length 


Lo 


2000 km 


T 


Temperature 


To 


6000 K 


P 


Density 


Po 


6.0 X 10-^ kg m-=^ 


P 


Pressure 


Po = poToKB/mn 


2.97 N m-2 


V 


Velocity 


Vo = {ToKBlmnf''' 


7.035 km s'^ 


B 


Magnetic field 


Bo = {P'oPoToKB/iriHy^^ 


19.3 G 


t 


Time 


to = Lo(ro«B/mjj)-V2 


284 s 
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Fig. 1. — Initial magnetic field configuration. The range of the computational box is from 
-1000 km to 1000 km in the x direction and from to 2000 km in the y direction. 
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Fig. 2. — Initial distributions of the temperature (left panel), the number density of hydrogen 
(middle panel) , and the electron number density (right panel) . The distributions in the quiet- 
Sun atmospheric model VALC (dashed lines) are also shown. 
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Fig. 3. — Temperature distributions (color scale), projected magnetic field (solid lines) and 
velocity field (vector arrows) at t = and t = 0.4 (upper panels) and t = 0.8 and t = 1.1 ta 
(lower panels). 
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Fig. 4. — Evolution of the magnetic reconnection rate \Vin\/VA- The rate reaches the maxi- 
mum at the time of 1.1 ta (120.8 s). 




Fig. 5. — One-dimensional T plot (solid line) along the y-axis at the time of the maximum re- 
connection rate (1.1 ta)- The initial temperature plot (dashed line), i.e. VALC temperature 
distribution, is also shown. 
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Fig. 6. — One-dimensional Vy plot (solid line) along the y-axis at the time of maximum 
reconnection rate (1.1 ta)- The height of the X-point is marked by a vertical line. The 
initial Vy plot (dashed line) is also shown. 
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Fig. 7. — Z-component of the current density (left panel) and temperature distributions 
(right panel) at the time of the maximum reconnection rate (1.1 ta). There is no need to 
plot the color bar for the left panel which based on dimensionless data, but in both panels 
the same color table is used, which means that the red side stands for higher current density 
and temperature, while the white side stands for the lower ones. 




Fig. 8. — Evolution of the magnetic reconnection rate for different initial background mag- 
netic fields (B) and anomalous resistivity (770). The height of the initial X-point is 1000 
km. 
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Fig. 9. — Temperature distributions for different initial background magnetic field (B) and 
anomalous resistivity {rjo) at the time of maximum rate of magnetic reconnection which have 
been shown in Figure IH1 The height of the initial X-point is 1000 km. 
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Fig. 10. — Temperature enhancement per unit time as a function of the background magnetic 
field B (left panel), the anomalous resistivity r^o (middle panel), and the number density of 
hydrogen atom uh (right panel). In the left panel, the initial reconnection height is 1000 
km and the anomalous resistivity is 17.7 Qm. In the middle panel, the initial reconnection 
height is 1000 km and the magnetic field is 25 G. In the right panel, the anomalous resistivity 
is 17.7 VLto. and the magnetic field is 25 G. 
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Fig. 11. — Magnetic reconnection rate (left panel), T (middle panel) and Vy (right panel) 
plot (solid line) at the time 0.24 ta- The initial conditions of T and Vy are also shown 
(dashed lines). The height of X-point is marked by the vertical line in the right panel and 
the temperature distribution of bright microflare (MF) is given by dash-dotted line in the 
middle panel. 
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Fig. 12. — Z-component of current density (left panel) and temperature distributions (right 
panel) at the time of maximum reconnection rate (0.24 r^). Both panels have used the same 
color table, which means that the red side stands for higher current density and temperature, 
the white side stands for the lower ones. 



